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We derive high-temperature series expansions for the free energy and the susceptibility of random- 
bond g-state Potts models on hypercubic lattices using a star-graph expansion technique. This 
method enables the exact calculation of quenched disorder averages for arbitrary uncorrelated cou- 
pling distributions. Moreover, we can keep the disorder strength p as well as the dimension d as 
symbolic parameters. By applying several series analysis techniques to the new series expansions, 
one can scan large regions of the (p, d) parameter space for any value of q. For the bond-diluted 
4-state Potts model in three dimensions, which exhibits a rather strong first-order phase transition 
in the undiluted case, we present results for the transition temperature and the effective critical 
exponent 7 as a function of p as obtained from the analysis of susceptibility series up to order 18. A 
comparison with recent Monte Carlo data (Chatelain et al, Phys. Rev. E64, 036120(2001)) shows 
signals for the softening to a second-order transition at finite disorder strength. 

PACS numbers: 

05.50.+q Lattice theory and statistics (Ising, Potts, etc.) 

64.60.Cn Order-disorder transformations; statistical mechanics of model systems 
64.60.Fr Equilibrium properties near critical points, critical exponents 



I. INTRODUCTION 

Systematic series expansions for statistical models de- 
fined on a lattice are a well-known method to study phase 
transitions and critical phenomena Q . They provide an 
useful complement to large-scale numerical simulations, 
in particular for quenched, disordered systems where the 
average over many different disorder realizations is nu- 
merically very time consuming and only some points in 
the vast parameter space of the systems can be sampled 
with realistic effort. 

Using high-temperature series expansions, on the other 
hand, one can obtain for many quantities results which 
are exact up to a certain order in the inverse tempera- 
ture. Here the infinite- volume limit can be taken without 
problems and the quenched disorder is treated exactly. 
Moreover, one can keep the disorder strength p as well 
as the dimension d as symbolic parameters and there- 
fore analyse much larger regions of the parameter space 
of disordered systems. To this end we developed further 
the method of "star-graph expansion" which allows us 
to take the disorder average on the level of individual 
graphs exactly and apply it to g-state Potts models with 
a bimodal quenched distribution of ferromagnetic cou- 
plings. Using optimized cluster algorithms for the exact 
calculation of spin-spin correlators on graphs with arbi- 
trary inhomogeneous couplings, we obtained series up to 
order 18 in the inverse temperature for the susceptibility 
and the free energy of bond-diluted Potts models in two, 
three and four dimensions. 
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Depending on the dimension d and the number of 
states q, pure Potts models show first- or second-order 
phase transitions. According to the Harris criterion Q 
one expects in the second-order case either the appear- 
ance of a new random fixed point (d = 2, q = 3, 4 and 
d = 3, q = 2) or logarithmic corrections to the pure fixed 
point (d = 2, q = 2). At first-order transitions, random- 
ness softens the transitions. For d = 2 even infinitesimal 
disorder induces a continuous transition ||, whereas for 
d — 3, q > 2 a tricritical point at a finite disorder strength 
is expected [Q. This softening to a second-order phase 
transition beyond a tricritical point at some finite disor- 
der strength has recently been verified in Monte Carlo 
(MC) simulations of the three-dimensional sife-diluted 
3-state H and frond-diluted 4-state || Potts model. 

The critical part of the series expansion methods lies in 
the extrapolation techniques which are used in order to 
obtain information on the critical singularity from a finite 
number of known coefficients of the high-temperature se- 
ries. One can question the use of these extrapolation 
techniques in disordered systems, where the complete 
singularity structure of the function may be very compli- 
cated, involving Griffiths-type singularities or logarith- 
mic corrections 

Anyhow, we are able to determine the transition tem- 
perature for the bond-diluted 4-state Potts model in 
three dimensions reliably up to the vicinity of the ge- 
ometrical percolation point and in good agreement with 
analytic estimates [|] and MC results Q. 

The critical exponent 7 extracted from our analysis ap- 
pears to be dependent on the disorder strength which is 
caused by crossover effects and the complicated singular- 
ity structure. Using sophisticated analysis methods, we 
find a range of the disorder strength where 7 c fr = 1, in- 
dicating critical behavior governed by a tricritical point. 

The rest of the paper is organized as follows. In Sect. II 
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we briefly recall the model. In Sect. Ill we describe the 



methods used for generating the series, and Sect. |TV| is de- 
voted to a representation of the analysis techniques used 
and their application to the study of the bond-diluted 
4-state Potts model in three dimensions. 



II. MODEL 



An articulation point is a vertex the deletion of which 
renders the graph disconnected. A graph without artic- 
ulation points is called "star graph" . 

A graph is bipartite if the vertices can be separated 
into red and black vertices so that no link connects two 
vertices of the same color. Equivalently, all closed paths 
in the graph consist of an even number of links. 



The g-state Potts model on the hypercubic lattice Z d , 
or more generally on any graph G with arbitrary coupling 
constants Jij assigned to the links (ij) of G, is defined 
by its partition function 



(1) 



where (3 = 1/kgT is the inverse temperature, Si — 
1, . . . , q and S(., .) is the Kronecker symbol. In our series 
expansion the combination 
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will be the relevant expansion parameter. In the sym- 
metric high-temperature phase, the susceptibility cor- 
responding to the coupling to an external field hi, 

hi ( qS ( Sl ^> — - J ( i s given for a graph with N spins by 
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Quenched disorder averages [. . . ] are taken over an un- 
corrected bimodal distribution of the form 

P{ J^) = (1 - p)5(J ij - Jo) + p5(Jij - RJ ), (4) 

which can include spin glasses (R = —l,p = 1/2), 
random-bond fcrromagnets (0 < R < 1) and bond di- 
lution (R = 0) as special cases. Other distributions can, 
in principle, also be considered with our method. 



III. SERIES GENERATION 

A. Basic notations from graph theory 

A graph of order E consists of E links connecting N 
vertices. We consider only connected, undirected graphs 
that are simple: no link starts and ends at the same ver- 
tex and two vertices are never connected by more than 
one link. Subgraphs are defined by the deletion of links. 
In this process, isolated vertices can be dropped. A graph 
of order E has 2 E (not necessarily non-isomorphic) sub- 
graphs since each link may be present or absent. These 
subgraphs may consist of several connected components 
and are called clusters. 



B. Star-graph expansion method 

There are two well-established methods |lj for the 
systematic generation of high-temperature series expan- 
sions, the linked cluster and the star-graph method. The 
longest known series (up to order (3 25 ) for classical spin 
models without disorder are produced by linked cluster 
expansions This technique allows one to obtain se- 
ries for observables (such as the second moment of the 
spin-spin correlation function) which have no star-graph 
expansion. Furthermore, it works with free embeddings 
of graphs into the lattice which can be counted orders of 
magnitude faster than the weak embedding numbers used 
by the star-graph technique. Nonetheless, this method 
has not yet been applied to problems with quenched dis- 
order. 

The star-graph method can be adopted to systems in- 
volving quenched disorder [jl0|, |ll| since it allows one 
to take the disorder average on the level of individual 
graphs. The basic idea is to assemble the value of some 
extensive thermodynamic quantity Fona large or even 
infinite graph from its values on subgraphs: Graphs con- 
stitute a partially ordered set under the "subgraph" rela- 
tion. Therefore, for every function F(G) defined on the 
set of graphs exists another function Wf(G) such that 
for all graphs G 



F(G) = J2 W F (g), 

gCG 

and this function can be calculated recursively via 
W F (G)=F(G)-^W F (g). 

9CG 

This gives for an infinite (e.g. hypercubic) lattice 
F{Z d ) =J2( G : ^ d )W F (G), 

G 



(5) 



(6) 



(7) 



where (G : Z d ) denotes the weak embedding number of 
the graph G in the given lattice structure [12] . 

The following observation makes this a useful method: 
Let G be a graph with an articulation vertex where two 
star subgraphs G\^ are glued together. Then W F (G) 
vanishes if 



F(G) — F(Gi) + F(G 2 ) 



(8) 
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TABLE I: Number of star graphs with E links and non- 
vanishing embedding numbers on Z d . 
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An observable F for which Eq. (^) is true on arbitrary 
graphs with articulation points allows a star-graph ex- 
pansion. All non-star graphs have zero weight Wf in the 
sum Eq. @). 

It is easy to see that the (properly normalized) free 
energy log Z has this property and it can be proved ficfl 
that the inverse susceptibility 1/x has it, too, even for 
arbitrary inhomogeneous couplings Jij . This restricts the 
sum in Eq. (Q) to a sum over star graphs. The linearity of 
Eqs. (||)-(7]) enables the calculation of quenched averages 
over the coupling distribution on the level of individual 
graphs. 

The resulting recipe for the susceptibility series is: 

• Graph generation and embedding number count- 
ing. 

• Calculation of Z{G) and the correlation matrix 
M nm (G) = Tr(q5(S n ,S m ) - l) e -/^«^» 

for all graphs as polynomials in E variables Vij. 

• Inversion of the Z polynomial as a series up to the 
desired order. 

• Averaging over quenched disorder, 
N nm (G) = [M nm /Z] p ^j^ , 

resulting in a matrix of polynomials in (p, v) . 

• Inversion of the matrix N nm and subgraph subtrac- 
tion, 

W X (G) = T, n , m (N- 1 ) nm - E gC G W x (g). 

• Collecting the results from all graphs, 
l/x = E G (G:Z d ) W X (G). 



C. Generation of star graphs and calculation of 
embedding numbers 

The most complicated part in every attempt to gener- 
ate lists of graphs by recursively adding nodes and edges 
to a smaller list is the isomorphism test, i.e., the decision 
whether two given adjacency lists or adjacency matrices 
describe the same graph modulo relabelling and reorder- 
ing of edges and nodes. We used the nauty package by 
McKay |l3[ ] which makes very fast isomorphism tests by 
calculating a canonical representation of the automor- 
phism group of the graphs. 

Since we are only interested in star graphs with non- 
vanishing weak embedding numbers in Z d , the following 
simple observations are helpful: 

• Only bipartite graphs occur since Z d is bipartite. 



12048(3) 
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FIG. 1: Two star graphs of order 17 and 19 and their weak 
embedding numbers up to 6 dimensions. 



A generic fc-dimensional embedding (i.e. one which 
really needs all k dimensions) contributes in d > k 
dimensions with degeneracy (f) . 



• A biconnected graph of odd order E = 2n + 1 has 
generic embeddings only up to dimension n since it 
must have at least two edges in each dimension. 

• The only biconnected graph of even order E = 2n 
which has generic embeddings of dimension n is the 
cycle of length 2n. All the other graphs will use at 
most 7i—l dimensions. 

For the embedding count we implemented a refined 
version of the backtracing algorithm by Martin |l2] . We 
did extensive tests to find the optimal algorithm for the 
"innermost" loop, the test for collisions in the embed- 
ding, and ended up using optimized hash tables. 

By this means, we classified for the first time all star 
graphs up to order 19 which can be embedded in hypercu- 
bic lattices (see Table |) and calculated their (weak) em- 
bedding numbers for (i-dimensional hypercubic lattices 
(up to order 17 for arbitrary d, order 18 and 19 for di- 
mensions < 4), see Fig. [j] for typical results. 



D. Cluster representation 

The partition function and the matrix of correlations 
M nm for each graph are calculated with arbitrary sym- 
bolic couplings Jij using the cluster representation 



Z = q N - E Y[(e^-l + q)Z, 



Z = q w Tr^Q [1 - + v lj qS(S t , Sj)] 



(9) 



(10) 
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Here the sum goes over all clusters C C G, E is the 
number of links (= order) of the graph G, e is the num- 
ber of links of the cluster and c the number of connected 
components of C. Z is normalized such that log Z has a 
star-graph expansion. This essentially reduces the par- 
tition sum from a sum over q N states to a sum over 2 E 
clusters. In the Ising case q = 2 another huge simplifica- 
tion takes place since only clusters where all vertices are 
of even degree contribute to the cluster sum. 

The 2 E clusters belonging to a graph are enumerated 
by Gray codes [jl4j such that two consecutive clusters in 
the sum (|ll]) differ by exactly one (added or deleted) link. 
Gray codes are a reordering of the binary representation 
of numbers such that the difference to the successor is in 
exactly one bit position. For example, for E = 4 the se- 
quence is 0000, 0001, 0011, 0010, 0110, 0111, 0101, 0100, 
1100, 1101, 1111, 1110, 1010, 1011, 1001, 1000 where ze- 
ros denote the deleted links. This allows to speed up the 
calculation considerably by re-using every term in the 
sum for the calculation of the next one. 

The calculation of the susceptibility involves the ma- 
trix of correlations M nm . The effect of inserting 
qS ( s ^S]) — 1 m t the trace of Eq. ( |l0| ) can easily be seen: 
we get one if the vertices n and m belong to the same 
connected component of the cluster and zero otherwise. 
Therefore, 



M nm oc ) _ q 




n (!-%•)), 



(12) 



where the sum is restricted to all clusters C nm C G in 
which the vertices n and m are connected. 

For the symbolic calculations we developed a C++ 
template library using an expanded degree-sparse rep- 
resentation of polynomials and series in many variables. 
The open source library GMP is used for the arbitrary- 
precision arithmetics. 

Our longest series, up to order 18, are obtained for the 
case of bond dilution where (Q) simplifies to 

P{J l3 ) = (1 - p)5(J ij - Jo) + p8{Jij), (13) 

since in this case the disorder average for a series is most 
easily done via 



„, n k 



«J*]p(j) = (i-p)X 1 



k„,n 1 + ...+n k 



(14) 



IV. 



SERIES ANALYSIS: TECHNIQUES AND 
RESULTS 



In the following we shall illustrate the analysis using 
the bond-diluted 4-state Potts model in three dimensions 
as our primary example. This model exhibits in the pure 
case a strong first-order transition which is expected to 
stay first order up to some finite disorder strength, before 
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FIG. 2: Ratio approximants for different dilutions p vs. 1/n. 



it gets softened to a second-order transition governed by 
a disordered fixed point. 

In the latter case we are interested in locating power- 
law divergences in the susceptibility series of the form 



X(v) = A(v c 



(15) 



For such a critical behavior many different series analysis 
techniques have been discussed in the literature which all 
have their merits and drawbacks |D| . 

To localize a first-order transition point, however, a 
high-temperature series alone is not sufficient since there 
the correlation length remains finite and no critical sin- 
gularity occurs. In analysing series by ratio, Pade or 
differential approximants, the approximant will provide 
an analytic continuation of the thermodynamic quanti- 
ties beyond the transition point into a metastable region 
on a pseudo-spinodal line with a singularity T* < T c and 
effective "critical exponents" at T*. 

Employing the techniques described above we obtained 
the high-temperature series expansions for the suscep- 
tibility up to order 18 with coefficients given as poly- 
nomials in the disorder strength p, as listed in the Ap- 
pendix for general dimensions d < 4. For such a series in 
two variables, the method of partial differential approx- 
imants |l6) should be well suited. Up to date, however, 
the only application of this method to a tricritical point 
p7[ used a test series of order 50 generated from an ex- 
actly solvable model. In our case, it was unable to give 
conclusive results. Therefore, we confined ourselves to 
the analysis of single-parameter series for selected values 
of p. 

The ratio method is the least sophisticated method of 
series analysis, but usually it is quite robust and gives 
a good first estimate of the series behavior. It assumes 
that the expected singularity of the form ([l5|) is the clos- 
est to the origin. Then the consecutive ratios of series 
coefficients behave asymptotically as 
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(16) 
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FIG. 3: Critical temperature for different dilution p as ob- 
tained from MC simulations Q and DLog-Pade series analy- 
ses. The inset shows the difference between the two estimates. 

Figure ^| shows these ratios for different values of p. In or- 
der to make them visually comparable, they are normal- 
ized by their respective critical couplings v c . For small 
p they show the typical oscillations related to the exis- 
tence of an antiferromagnetic singularity at — v c . Near 
the percolation threshold at p — 0.751 188 (l8) (where 
T c goes to 0) the series is clearly ill-behaved, related to 
the exp(l/T) singularity expected there. Besides that, 
we observe that the slope (oc 7 — 1) is increasing with p, 
changing from 7 < 1 to 7 > 1 around p = 0.5. 

The widely used DLog-Pade method consists in calcu- 
lating Pade approximants to the logarithmic derivative 
of x(v). The smallest real pole of the approximant is an 
estimation of v c and its residue gives 7. Figure ^ com- 
pares the critical temperature, estimated from an average 
of 25-30 Pade approximants for each value of p p3] , with 
the results of recent MC simulations || . For small p, in 
the first-order region, the series underestimates the criti- 
cal temperature. As explained above, this is an estimate 
not of T c but of T*. Between p = 0.3 and p = 0.5, the 
estimates confirm, within errors, the MC results, indi- 
cating that now both methods see the same second-order 
transition. Beyond p = 0.5, the scatter of different Pade 
approximants increases rapidly, related to the crossover 
to the percolation point. 

The situation is more complicated with respect to the 
critical exponent 7. A DLog-Pade analysis gives incon- 
clusive results due to a large scattering between different 
Pade approximants, as shown in Fig. ||. One possible 
reason for this failure is the existence of confluent singu- 
larities: The dots in Eq. ( p"5| ) indicate correction terms 
which can be parametrized as follows: 

X (v) = A{v c - «)-T[l + A x {v c - v) Al + A 2 (v c - v) A2 +... 

(17) 

where Aj are the confluent correction exponents. We 
used different more sophisticated analysis methods, such 
as inhomogeneous differential approximants Jig] and the 
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FIG. 4: Scattering of different Pade approximants at a di- 
lution p — 0.4: critical exponent 7 against critical coupling 

Vc- 

methods Ml and M2 p9| , especially tailored to deal with 
confluent singularities. In the case at hand, the Baker- 
Hunter method |2(J appeared to be quite successful, giv- 
ing consistent results at larger dilutions p > 0.35 where 
the DLog-Pade analysis failed. Assume the function un- 
der investigation has confluent singularities 

N , v-Aj 

F(z)=X> 1-- = ^a n z n . (18) 

j=l ^ c ' n=0 

This can be transformed into an auxiliary function g(t) 
which is meromorphic and therefore suitable for Pade 
approximation. After the substitution z — z c (l — e ) 
we expand F(z(t)) = ^2 n c n t n and construct the new 
series 

N , 

g (t) = Y / n\c n t n = J2 T zj-t ( 19 ) 

n=0 i=l 1 

We see that Pade approximants to g(t) have poles at 
t = l/Xi, with residues at the poles of — AijXi. This 
method is applied by plotting these poles and residues 
for different Pade approximants to g(t) as functions of 
z c . The optimal set of values for the parameters is deter- 
mined visually from the best clustering of different Pade 
approximants, as demonstrated in Fig. ^. 

Using this method, our results for the critical exponent 
7 are plotted in Fig. |^. They show an effective exponent 
monotonically increasing with p but reaching a plateau at 
7=1 for dilutions between p = 0.42 and p — 0.46. The 
following sharp increase is to be interpreted as due to the 
crossover to the percolation fixed point p c = 0.751 188, 
T c = 0, where a x ~ exp(l/T) behavior is expected. 
; It is well known (see, e.g., Ref. ^T|) that series analysis 
in crossover situations is extremely difficult. If the pa- 
rameter p interpolates between regions governed by dif- 
ferent fixed points, the exponent obtained from a finite 
number of terms of a series expansion must cross some- 
how between its universal values, and does this usually 
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FIG. 5: Values for the critical exponent 7 and amplitude A at 
p — 0.4 as function of trial v c estimates from the Baker-Hunter 
analysis. From the clustering of different Pade approximants 
in both pictures we estimate v c = 0.3217, 7 = 0.966, and 
A = 1.21. 
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FIG. 6: Effective critical exponent 7 as function of the dilu- 
tion p from Baker-Hunter analyses. 

one can obtain this average exactly. Since the relevant 
parameters (degree of disorder p, spatial dimension d, 
number of states q, etc.) can be kept as symbolic vari- 
ables, the number of potential applications is very large. 

Here we presented an analysis of the three-dimensional 
bond-diluted 4-state Potts model. The phase diagram 
confirms recent Monte Carlo results and, by comparing 
with the numerical data, we also see signals for the onset 
of a second-order transition at a finite disorder strength. 
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quite slowly. Therefore it does not come as a surprise 
that the MC simulations quoted above see the onset of 
a second order phase transition already for smaller val- 
ues of the disorder strength p. The mere existence of 
a plateau in 7 e ff(p), however, is an indication that here 
truly critical behavior is seen. It is governed by a fixed 
point for which we obtain 7 = 1.00(3). Here, as always 
in series analyses, the error estimates the scattering of 
different approximants. 



V. CONCLUSIONS 



We have implemented a comprehensive toolbox for 
generating and enumerating star graphs as required for 
high-temperature series expansions of quenched, disor- 
dered systems. Monte Carlo simulations of systems with 
quenched disorder require an enormous amount of com- 
puting time because many realizations have to be sim- 
ulated for the quenched average. For this reason it is 
hardly possible to scan a whole parameter range. Using 
high-temperature series expansions, on the other hand, 
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APPENDIX A 

As an example, we publish here the inverse susceptibility for the bond-diluted 4-state Potts model in d < 4 
dimensions up to order v 18 (P = 1 — p): 

X _1 (P,v,d) = l-2Pvd + 2P 2 v 2 d-2P' A v :i d+ [2P 4 d~ 16P 4 (£)] v A + [-2 P 5 d + (24 P 4 + 72 P 5 ) (£)] t> 5 
+ [2P 6 d+(24P 4 -96P 5 - 248P 6 ) $ - 768P 6 (3)] v 6 

+ [-2P 7 d + (24 P 4 -96P 5 + 264P 6 + 640P 7 ) (J + (576 P 6 + 3264 P 7 ) Q] v 7 

+ [2P 8 d + (72 P 4 - 96 P 5 +264P 6 - 216 P 7 - 1384 P 8 ) (£) + (576 P 6 - 144 P 7 - 22704 P 8 ) (f) - 62208 P 8 (f)] v 8 

+ [-2 P u d + (-72 P 4 - 288 P 5 + 264 P 6 - 312 P 7 - 1416 P 8 + 1888 P 9 ) (£) 
+ (576 P 6 - 720 P 7 + 720 P 8 + 66944 P 9 ) (£) + (31104 P 8 + 221312 P 9 ) (*)] v 9 

+ [2 P 1 " d + (-72 P 4 + 288 P 5 + 648 P 6 + 384 P 7 - 144 P 8 + 9336 P 9 - 296 P 10 ) $ 

+ (576 P 6 + 3456 P 7 + 19296 P 8 + 75456 P 9 - 387168 P 10 ) (f) + (31104 P 8 + 109440 P 9 - 4000512 P 10 ) (f)] v 10 

+ [-2P"d + (-72 P 4 + 288 P 5 - 504 P 6 - 1008 P 7 - 3024 P 8 + 3144 P 9 - 33336 P 10 - 9616 P 11 ) (J 
+ (576 P 6 - 1440 P 7 - 8352 P 8 - 31248 P 9 - 309744 P 10 + 781824 P 11 ) Q 
+ (31104 P 8 + 11520 P 9 + 635520 P 10 + 10415872 P 11 ) (^)] u 11 

+ [2 P 12 d + (-216 P 4 + 288 P 5 - 216 P 6 - 4272 P 7 + 240 P 8 + 11856 P 9 - 4968 P 10 + 81744 P 11 + 37320 P 12 ) (J 
+ (2880 P 6 - 31392 P 7 + 14112 P 8 + 169200 P 9 + 489024 P 10 + 2692800 P 11 - 5811664 P 12 ) Q 
+ (31104 P 8 + 273024 P 9 + 3204864 P 10 + 16037760 P 11 - 179275648 P 12 ) (f)] v 12 

+ [-2 P 13 d + (216 P 4 + 864 P 5 - 792 P 6 - 3912 P 7 + 29736 P 8 + 5952 P 9 - 20736 P 10 + 23088 P 11 
-144624 P 12 - 96160 P 13 ) (J 

+ (-1728 P 6 - 31536 P 7 + 162288 P s + 15408 P 9 - 223344 P 10 - 113760 P 11 - 8412192 P 12 + 5990784 P 13 ) Q 
+ (31104 P s +67968 P 9 + 1022976 P 10 - 693504 P 11 - 16255872 P 12 + 304010112 P 13 ) (f)] v rs 

+ [2 P 14 d + (216 P 4 - 864 P 5 - 1944 P 6 - 8616 P 7 + 18360 P 8 - 99600 P 9 - 65544 P 10 + 33936 P 11 
-86952 P 12 + 73704 P 13 + 169400 P 14 ) (£) 

+ (-1728 P 6 - 55152 P 7 + 67248 P s - 894240 P 9 - 918000 P 10 + 2799648 P 11 + 8589744 P 12 
+58983984 P 13 - 98045424 P 14 ) (3) 

+ (31104 P s - 1057536 P 9 + 17280 P 10 + 24870528 P 11 + 179980416 P 12 + 1095494784 P 13 - 7487817088 P 14 ) (f)] v 14 

+ [-2P 15 d+ (216 P 4 - 864 P 5 + 1512 P 6 - 4536 P 7 + 44568 P 8 - 55200 P 9 + 168480 P 10 + 363072 P 11 - 11832 P 12 
+530040 P 13 + 501600 P 14 - 145632 P 15 ) (^) 

+ (-1728 P 6 - 41040 P 7 + 209232 P s - 609984 P 9 + 1319328 P 10 + 7874208 P 11 + 5670048 P 12 + 24319296 P 13 
-141840288 P 14 - 14817536 P 15 ) (3) 

+ (31104 P 8 - 740736 P 9 + 131328 P 10 + 22334976 P 11 + 66366720 P 12 + 188319744 P 13 
-1467511296 P 14 + 5362518016 P 15 ) (f)] -y 15 

+ [2P 16 d+ (648 P 4 - 864 P 5 + 1512 P 6 + 24336 P 7 + 38496 P 8 - 148008 P 9 + 1656 P 10 + 117024 P 11 - 1325376 P 12 
+7200 P 13 - 1644000 P 14 - 2926176 P 15 - 373984 P 16 ) (^) 

+ (-1728 P 6 + 163296 P 7 + 105984 P 8 - 2305728 P 9 - 1857888 P 10 - 6620544 P 11 - 48148992 P 12 + 18163728 P 13 
+118520640 P 14 + 1144225008 P 15 - 1918717248 P 16 ) (3) 

+ (217728 P 8 - 3438720 P 9 - 7119360 P 10 - 64137600 P 11 - 149601024 P 12 + 1152714240 P 13 + 8368094208 P 14 
+58294742400 P 15 - 317165909504 P 16 ) (f)] « 16 

+ [-2 P 17 d + (-648 P 4 - 2592 P 5 + 1512 P 6 + 15408 P 7 - 161712 P 8 - 123384 P 9 + 66792 P 10 + 39264 P 11 
-1976760 P 12 + 3413424 P 13 + 848256 P 14 + 4241568 P 15 + 8541960 P 16 + 2398960 P 17 ) (3) 

+ (-1728 P 6 + 109728 P 7 - 984096 P s - 1584432 P 9 + 413424 P 10 - 11287296 P 11 - 32069376 P 12 + 206240976 P 13 
+380730960 P 14 + 1087235856 P 15 - 1859056704 P 16 - 2643006384 P 17 ) (3) 

+ (-93312 P s - 2522880 P 9 - 4468608 P 1() - 63930240 P 11 - 127255680 P 12 + 1693209600 P 13 + 6161021568 P 14 
+23385824256 P 15 - 50368269312 P 16 - 105383991680 P 17 ) (f)] v 17 
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+ [2 P 18 d + (-648 P 4 + 2592 P 5 + 4104 P 6 + 83520 P 7 - 82080 P 8 + 465984 P 9 + 586080 P 10 + 605064 P 11 - 166248 P 
+8121312 P 13 - 4714536 P 14 - 2886168 P 15 - 3604536 P 16 - 20651832 P 17 - 7297424 P 18 ) (J 

+ (-8640 P 6 + 518400 P 7 - 492480 P 8 + 3752592 P 9 + 12513744 P 10 - 23522256 P 11 - 66704640 P 12 - 53106912 P 13 
-884626272 P 14 - 122168448 P 15 + 2205877392 P ie + 22700601216 P 17 - 42014019168 P 18 ) (3) 
+ (-93312 P 8 + 657024 P 9 + 23180544 P 10 - 169350912 P 11 - 762268032 P 12 - 3977024256 P 13 - 10126195200 P 14 
+47139877632 P 15 + 379559824128 P 16 + 2866361546496 P 17 + 3747410465664 P 18 ) (^] w 18 
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